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AERONAUTIC  SYMBOLS 
1.  FUNDAMENTAL  AND  DERIVED  UNITS 


i 


w 

e 

m 

1 

S 

O 

h 

c 

A 

V 

2 
L 
D 
Do 
Di 

0 


Symbol 

Metric 

English 

Unit 

Abbrevia¬ 

tion 

Unit 

Abbrevia¬ 

tion 

1 

t 

F 

m 

8 

kg 

foot  (or  mile) _ 

ft  (or  mi) 
sac  (or  hr) 
lb 

Time _ 

Force - 

second _ 

weight  of  1  kilogram _ 

second  (or  hour) _ 

weight  of  1  pound _ 

P 

V 

horsepower  (metric) _ 

horsepower _ 

lip. 

mph 

fps 

Speed - 

fldlomieters  per  hour _ 

\metera  per  second - 

CD 

miles  per  hour _ 

feet  per  second _ 

2.  GENERAL  SYMBOLS 


Weight==mi; 

Standard  acceleration  of  ^avity= 9.80665  la/s’ 
or  32.1740  ft/secV  :  , 

Masa=^ 

Moment  of  inertia=ml;*.  (Indicate  axis  of 
radius  of  gyration  k  by  proper  subscript.) 
Coefficient  of  viscosity 


r  Kinematic  viscosity 

p  Density  (mass  per  unit  volume) 

Standard  density  of  dry  air,  0.12497  kg-m“^-s*  at  16®  C 
and  760  mm;  or  0.002378  Ib-ft"*  sec* 

Specific  weight  of  “standard”  air,  1,2265  kg/m*  or 
0.07661  Ib/cu  ft 


8.  AERODYNAMIC  SYMBOLS 


Area 

Area  of  wing  - 

Gap 

Span 

Chord 

Aspect  ratio, 

True  air  speed 
Dynamic  pressure,  |pF* 

Lift,  absolute  coefficient  Cl=^ 

Drag,  absolute  coefficient  Cd—^ 


Do 


Profile  drag,  absolute  coefficient  Cc(,=^ 


Induced  drag,  absolute  coefficient 

^  go 


Parasite  drsLg,  absolute  coefficient 

<7 

Cross-wind  force,  absolute  coefficient 


iw 

it' 

Q 

a 

R 


a 

e 

“0 

o< 

Oo 
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Angle  of  setting  of  wings  (relative  to  thrust  line) 
Angle  of  stabilizer  setting  (relative  to  thrust 
line)  '  ;  I 

Resultant  moment 
Resultant  angular  velocity 


Re3molds  number,  p—  where  1  is  a  linear  dimen¬ 
sion  (e.g.,  for  an  airfoil  of  1 .0  ft  chord,  100  mph, 
standard  pressure  at  15®  C,  the  corresponding 
Reynolds  number  is  935,400;  or  for  an  airfoil 
of  1.0  m  chord,  100  mps,  the  corresponding 
Reynolds  number  is  6,865,000) 

Angle  of  attack 
Angle  of  downwash 
Angle  of  attack,  infinite  aspect  ratio 
Angle  of  attack,  induced 
Angle  of  attack,  absolute  (measured  from  zero- 
lift  position) 

Fhght-path  angle 
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ISOLATED  AND  CASCADE  AIRFOILS  WITH  PRESCRIBED  VELOCITY  DISTRIBUTION 


By  Arthir  \V.  Goldstein  and  Meyer  .Ierison 


SUMMARY 

-l/i  exact  solution  of  the.  problem  of  desie/ninij  an  airfoil  with 
a  prescribeil  relocity  ditttribution  on  the  suction  surface  in  a 
ijicen  uniform  flow  of  an  incompressible  perfect  fluid  is  obtained 
by  replaciny  the  boundary  of  the  airfoil  by  cortices.  By  this 
decice,  a  method  of  solution  is  developed  that  is  applicable  both 
to  isolated  airfoils  and  to  airfoils  in  cascade.  The  conformal 
transformation  of  the  designed  airfoil  into  a  circle  can  then  be 
obtained  and  the  velocity  distribution  at  any  angle  of  attack 
computed.  Numerical  illustrations  of  the  method  are  given  for 
the  airfoil  in  cascade.  ^  7"^^^ — 

INTRODUCTION 

The  problem  of  increasing  the  output  per  stage  in  axial- 
flow  compressors  and  turbines  involves  the  use  of  high- 
solidity  (closely  spaced  blades)  stages  of  highly  cambered 
blades.  In  addition,  the  velocity  distribution  must  be  care¬ 
fully  selected  as  a  function  of  arc  length  along  the  airfoil 
(blade  section)  boundary  in  orrler  to  avoid  flow  separation  or 
c.xcessivcly  high  local  velocities. 

Several  methods  are  available  for  obtaining  an  airfoil  with 
a  prescribed  velocity  distribution.  The  methods  that  lead  to 
theoretically  exact  results  are  based  on  conformal-mapping 
tlicory.  (See  references  1  and  2.)  In  reference  Mutterpcrl 
extends  the  method  of  conformal  mapping  to  solve  the 
problem  of  computing  a  cascade  of  airfoils  with  prescribed 
velocity  distribution  but.  for  cascades  with  closely  spaced  or 
highly  cambered  airfoils,  this  procedure  becomes  very  cum¬ 
bersome.  Approximate  solutions  have  been  obtained  by 
placing  singularities  such  as  vortices,  sources,  and  sinks  in  a 
uniform  stream.  The  shape  of  sections  of  airfoils  in  cascade 
can  also  be  computed  by  distributing  such  singularities 
licriodically  throughout  the  region  of  the  cascade,  as  ilcscribed 
by  Ackcrct  (reference  4). 

Because  these  vortex  methods  arc  not  e.xact.  a  method 
witli  the  vortices  on  the  boundaries  of  the  cascade  airfoils 
was  developed.  This  method  gives  a  theoretically  exact  solu¬ 
tion  without  the  computation  difficulties  encountered  in 
conformal-mapping  methods  for  highly  cambered  airfoils  or 
closely  spaced  cascades.  Furthermore,  for  the  same  accuracy 
in  computing  the  airfoil  shape,  this  vortex  method  requires 
tlie  computation  of  fewer  points  than  the  mcthoil  of  conformal 
mapping  because  these  points  may  be  arbitrarily  placed  on 
the  airfoil.  The  method  may  be  applied  to  isolated  airfoils 
and  to  airfoils  in  cascade.  For  the  cascade,  the  inflow  and 
ilischarge  velocities  and  a  velocity  distribution  on  the  surface 
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of  an  airfoil  are  given  and  the  shape  of  the  airfoil  is  deter¬ 
mined.  In  some  cases,  the  spacing  of  the  blades  is  pre¬ 
assigned,  which  places  a  condition  on  the  assumed  velocity 
distribution.  Once  the  airfoil  shape  has  been  evolved,  the 
velocity  distribution  may  be  computed  for  any  angle  of 
attack  by  the  method  described  in  appendix  A.  The  method 
of  this  paper  was  developed  at  the  MACA  Cleveland  labora¬ 
tory  during  1946. 

THEORY 

OUTLINE  OF  METHOD 

In  reference  5,  it  is  demonstrated  that  the  two-dimensional 
potential  flow  about  a  body  in  a  uniform  stream  can  be 
represented  by  substituting  for  the  body  a  sheet  of  vortices 
along  its  boundary.  The  vortex  strength  per  arc  length  at 
any  point  is  equal  to  the  magnitude  of  the  velocity  at  that 
point.  A  proof  of  this  relation  for  the  case  of  the  cascade  is 
given  in  appendix  B.  The  problem  of  finding  a  shape  with 
a  prescribed  velocity  distribution  when  placed  in  a  stream 
can  then  be  stated:  Given  a  vortex  distribution,  to  find  a 
contour  which  satisfies  the  condition  that  it  will  be  a  stream¬ 
line  in  the  flow  field  induced  by  the  uniform  flow  and  the 
vortices  distributed  on  the  contour. 

The  procedure  of  finding  the  shape  begins  by  choosing  an 
approximate  shape  and  distributing  the  vortices  on  it.  The 
stream  function  of  the  flow  induced  by  the  vortices  and  the 
uniform  stream  is  computed  at  points  on  the  boundary  of 
the  assumed  shape.  If  this  stream  function  is  constant,  the 
assumed  shape  is  correct.  Variations  of  the  stream  function 
are  a  measure  of  the  deviation  of  the  assumed  shape  from 
the  correct  one.  These  variations  are  used  to  distort  the 
original  shape  into  a  new  shape  whose  stream  function  is 
more  nearly  constant.  The  process  is  repeated  until  the 
variations  become  negligible.  In  the  iirocess  of  shape  adjust¬ 
ment,  the  velocity  is  altered  on  the  pressure  surface. 

DERIVATION  OP  EQUATIONS  FOR  THE  STREAM  FUNCTION 

Isolated  airfoil. — The  complex  or  reflected  velocity 
w'(z)  (which  is  the  derivative  of  the  complex  potential  func¬ 
tion  w{z))  induced  at  the  point  z—xf-iy  by  a  vortex  of 
strength  k  located  at  2o  =  J.a-f-'?/o  is 


(A  summary  of  the  principal  symbols  used  in  this  report  is 
given  in  appendix  C.) 

1 


2 


IfEPOKT  XO.  S(>!) - NATIONAL  ADNTSORV  COMMITTKK  FOH  AERONAUTICS 


The  complex  velocity  )r'{:)  induced  hv  ii  uniform  slrenm 
with  comiilex  velocity  ir,,'  and  a  dislrihutioii  of  v'ortex 
strength  ])er  unit  lens;:th  7(r„)  alontr  a  curve  with  coordinates 

-  fli 

'  Cy(z,,)  <lx,. 


where  f/x„  is  the  element  of  arc  hmetth  along  the  curve.  The 
complex  iiotential  at  the  jioint  ;  is  the  integral  of  w'{z)  with 
respect  to  2,  namely 

=  +  I  7(2„)  log  (s-2„)  (/.%  (2) 

From  reference  1  (notation  modified), 

7(2„)  (/6'„=w'(2»)  dz^dw{z„)=dip{z„)  +  ldii{z„) 

where 

if  velocity  potential.  /?[«,'(2)] 

\p  stream  function,  I[w{z)\ 

Wien  ecpiation  (2)  is  applied  to  obtain  the  complex 
Iiotential  function  at  any  point  z  in  the  flow  field,  the  in¬ 
tegration  must  bo  carried  otit  along  the  hotindarv  of  the 
body.  Because  this  curve  is  a  streamline.  dip=0  and,  there¬ 
fore,  equation  (2)  heeomes 

U)(2)=2Wu'-f-2“-  J  log(2-2„)  (2a) 

The  imaginary  part  of  etpiation  (2a)  is  the  stream  function 
at  the  point  z, 


4>=-rV,+;/V,-,y^  j  log 


:  ^,  '{x-s„y+{y—y„yd,p{z„)  (3) 


whei'e 

V„  i/-component  of  tmiform  stream  velocity  F 
T’l  j-component  of  tmiform  stream  velocity  V 

It  is  convenient  to  use  the  arc  length  along  the  airfoil  as 
a  parameter.  If  (x.y)  is  a  point  on  the  airfoil  boundary,  then 
6'  will  denote  the  arc  length  there;  similarly,  will  denote 
the  arc  length  at  (j'„,  y„).  The  vortex  at  .v»  on  the  airfoil 


inllmmces  the  stream  fimi'tion  at  the  point  s  on  the  ai 
The  stream  function  induced  at  (./■.//)  by  a  vortex  of 
strength  at  (.r„.  //„)  is 


./'ll-, 


log 


A  plot  in  the  {.r.y)  idane  of  curves  for  a  constant  /', 
consists  of  concentric  circles  with  center  at  ix„,  y„). 

'riio  velocity  at  the  point  on  the  airfoil  is  the  direci 
derivative  p'(s„)  of  the  potential  along  the  strean 
If  the  velocity  along  the  airfoil  has  been  siiecified  an 
airfoil  shape  has  heim  assumed,  the  resultant  stream 
tion  along  the  houndary  of  the  airfoil  can  be  approxin 
by  using  the  aiiproximate  shape  in  evaluating  the  inleg 

'/'(■'>')  =  4',,  ('s)  —  jji  (■'■ ,  ■'•<,)  ‘fi'  (•'’■») 

where 

^„(s)  stream  function  at  (ji.,?/)dtie  to  uniform  stream,  — /F„- 
I  total  arc  length  of  airfoil 

All  variables  arc  expressed  in  term?  of  the  arc-length  i 
meters  s  and  .v„.  The  inti'gral  in  eqtiation  (5)  can  h(‘  e 
ated  either  numerically  or  graphically  over  the  entire  i 
of  integratioir  except  in  the  region  where  a  (s.<!— .<„)  is  s 
for  in  this  region  /i(.s',.s'„)  heeomes  infinite.  This  portic 
the  integral  can  be  evahiated  by  approximating  the  a 
hotindarv  hy  a  line  segment.  Then, 

«■»)  log  {s-s^y 

'Phe  ])r('seribed  velocity  can  he  given  in  this  ri'gion,  v 
may  be  defined  by  .s— « Ss„S6'-|-rt,  by  a  Taylor’s  seri, 
a  function  of  .<!„  about  the  point  s. 

v''(.sJ  =  v''(-)  +  V’''(''')(^'«-v)  +  ^^  (.v„-v)-+  •  •  ■ 

where  the  primes  indi<'ate  derivatives  with  respect 
The  integral  is  then 


Xj+o  P'+n  ] 

^ /i(s, ■'••„) ^  log  (.v„  — .s-)‘  [i,5'('V)  +  ',!’”('dl'-'''»— '-O-r  ■  ■  •] 

=  “[^<i,p'(.s-)  (log  «— l)-t---’“p—  (log  ■  -"I 

In  most  cases,  only  the  first  term  need  lie  used  in  equation  (ti).  The  same  type  of  ap|)roximation  can  he  used  to  evai 
a  portion  of  the  integral  if  the  opposite  side  of  the  airfoil  comes  in  the  neighhorhood  of  the  point  ix,y). 

A  more  general  e<iuation  applicable  to  a  segment  that  does  not  pass  through  s'  is: 


log  [(J'-J-„)--i-(y-i/„)-]¥’'(s„)  d(-v„)=^*^  '!«’'(/>)  j^c  log  {td+c-)-b  log  (/r'  +  i-) -2ic-6) -f  2/t  tau-‘  f 


[(/d-|-c-)  log  (ld-\-r)  —  \lryb'-)  log  (lr-\-b'-)  —  (c-  —  b-')]  + 


log  {h-+y)-b^  log  (y-P)  + 

2A^(c-6)_24Man-‘^;J^]+-  •  ■  ] 
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/(  is  tile  |)erpeu(liculai-  distiuiee  from  x  lo  (lie  seirnienl. 
.v„  =  /)  locates  tlio  foot  of  tlic  seirmeiit.  (/)+6f  and  f/i+e) 
are  the  liinils  of  tlie  inteirratioii  of  x,„  and  approximately, 


p  —  b  >—b^"'{p  —  li) 

I'Apialion  (lia)  may  lie  us(‘d  wluai  the  line  segment  is  nol 
of  equal  leiiftths  on  either  side  of  the  ix'rpendicidar  fool  or 
when  i^i'fx)  or  its  derivatives  are  diseontimioiis  at  either 
(li  +  b)  or  (/i  +  r).  If  «=(•=  — fi  and  /i=0,  e(|ualion  (lia) 
reduces  to  ('(piiition  (li).  'Phe  size  of  a,  b.  or  e  is  deter¬ 
mined  hy  the  requirements  that  the  septment  in  (piestion 
he  nearly  straight  (the  approximation  is  of  the  second  decree) 
and  that  <(>'{>>„)  he  accurately  represented  hy  a  'Paylor’s 
si'i'ies  expansion  of  few  terms. 


Exit-velocity  diagrarrt  Complete-velocity  diagram 


Entrance  -velocity  diagram  loz'  tom’  toi' 

Ref /ected- velocity  diagram 

Figure  1.— Nolnlion  for  c:w:ulc»  flow. 

Airfoils  in  cascade. — The  expression  for  the  complex  po¬ 
tential  for  the  flow  about  a  cascade  of  airfoils  is  derived  in 
appendix  B.  The  notation  is  defined  in  figure  1.  The 
eipiation  that  corresponds  to  equation  (2a)  for  isolated  air¬ 
foils  is  for  a  cascade  of  airfoils 

w(2)  =  I’log  l^sin  I  (.—  ;„)  J  r/,p(g„)  (7) 

where 

«'„/  mean  stream  velocity,  which  is  one-half  the  sum  of 
complex  (reflected)  velocities  upstream  and  down¬ 
stream  of  cascade.  — iV„ 

S  distance  betw'een  successive  airfoils  in  cascade 

The  mean  velocity  w„'  corresponds  to  the  uniform  velocity 

u'u  of  the  isolated  airfoil  flow. 

d'hc  term  zw^'  is  the  complex  potential  function  resulting 
from  the  mean  flow.  In  the  integral,  the  element  dip  indi¬ 


cates  the  vortex-element  strength  and  log  [sin  (iriS)  (---,)] 
represents  the  complex  potential  at  the  point  ;  caused  hy  an 
infinite  row  of  unit  vortici-s  at  .:„±mS'  where  »=(),  1, 
2,  .  .  ..  'Phe  imaginary  pari,  of  e(|uation  (71  is  the  stream 
function. 

'/'(■'')  =1/- —  (  fi(p.x.,)  <h(n^)  (8) 

where 

./2(x.x„)  =  log  j^siir’  J,  +siidi-  (//  — ,'/„)  J 

is  expressed  in  are-huigth  parameters  and  i/'m(.s’)  is  the  stream 
function  at  (r.y)  induced  hy  a  mean  stream  whose  complex 
velocity  is  w,„' ;  that  is, 

'Pm  =  V,j  -Tl/Vi 

'Phe  values  of  (.r-x„)/S  and  (y-i/of/S  for  various  values  of 
are  given  in  table  1.  A  plot  of  x-x,,  and  //-;/„  for  constant 
values  of  /■,{z,z„)  is  shown  in  figure  2.  'Phese  curves  may  bo 


y-Vol  (yVo) 


.3 


Figure  2.— Plot  of  curves  for  constant  fi(r— /a.  v-t/o). 

interpreted  as  the.  streamlines  of  the  flow  induced  by  an 
infinite,  row  of  vortices  of  unit  strength  located  at  the  points 
(t<,±hS.  i/„),  where  n=0,  1,  2,  .  .  ..  In  the  region  of  a 

vortex,  the  streamlines  are  nearly  circles;  that  is,  the  flow  is 
nearly  that  induced  by  an  isolated  vortex.  At  a  distance 
from  the  vortex  row,  the  streamlines  are  parallel  lines,  as  in 
the  flow  pattern  induced  by  a  continuous  uniform  distribu¬ 
tion  of  vorticity  along  a  straight  lino  instead  of  a  row  of 
discrete  vortices.  'Phe  velocities  on  the  two  sides  of  such  a 
vortex  line  are  of  eciual  magnitude  but  opposite  in  direction. 
This  behavior  of /o  for  largo  \y~>ji,\IS  and  also  for  small 

— -'i-  can  be  described  as  follows;  When  both 

{x.~x„)lS  and  {>j—yo)IS  are  small, 

Ms,  s,)  log  g-i  [(x-x^y  +  iy-ya)-]  (9) 

which  differs  from  /i(3,  z„)  only  by  a  constant.  For  largo 
values  of  ly—yai/S,  irrespective  of  (x—XojjS  and  a  constant 
term, 

(10) 

which  is  the  stream  function  of  a  uniform  stream  parallel 
to  the  x-a.xis. 
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Equation  tSt  can  he  used  for  computinn;  the  stream  function  aloii"  tlie  houndary  of  an  airfoil  in  cascade  j 
c(|uation  (o)  is  used  for  the  isolated  airfoil.  The  intefiral  over  the  raiifie  in  the,  neiirhhorhood  of  the' point  .v  is  obtaii 
using  equation  (9)  for/2  So).  The  result,  derived  in  the  same  maimer  as  equation  (6),  is 


The  more  general  e([uation  (ha)  is  modified  for  cascades  hy 
multiplying  the  argument  of  all  logarithms  hy  the  faetor 
r/S-. 

ADJUSTMENT  OF  SHAPE 

If  the  Stream  function  for  the  assumed  airfoil  has  been 
computed  and  has  been  found  to  vary,  the  shape  must  then 
he  adjusted  to  give  a  more  tiearly  constant  stream  function. 
The  shape  changes  are  tnade  by  rotation  of  the  body  plus 
displacement  of  the  individual  points  normal  to  the  mean 
stream.  The  rotation  is  used  to  place  the  front  stagnation 
point  properly. 

Rotation  of  the  airfoil. — In  the  formula  for  computing  the 
stream  function  of  an  isolated  airfoil,  the  contribution  of  a 
vortex  element  at  {Xo,  ]jo)  to  the  stream  function  of  a  point  at 
(/,  y)  is  dependent  merely  on  the  distance  between  the  two 
points.  Consequently,  if  the  entire  airfoil  is  rotated,  the 
effect  of  the  boundary  vortices  on  the  stream  function  at  any 
point  on  the  airfoil  houndary  will  not  change.  The  effect  of 
the  blade  rotation  on  the  stream  function  along  the  boundary 
is  therefore  determined  by  the  change  in  relative  position  of 
the  points  in  the  uniform  stream.  The  first  adjustment  in 
shape  is  a  rigid  rotation  of  the  airfoil  in  order  to  obtain  a 
more  nearly  constant  stream  function  along  its  boundary. 

If  the  airfoil  is  rotated  through  an  angle  d,  the  stream 
function  is  so  changed  that  xp  (x)  is  a  function  of  /3  and  s  and 


may  be  written  xp  (.v,  d).  When  ji  —  O.  xp  (.v,  0)  is  the  01 
stream  function  before  rotation.  After  rotation  tic 
stream  function  p  (.s,  d)  may  he  expanded  in  a  Taylor’s 
about  the  point  0  =  0, 

Pis,0)  =  P(s,O)  +  0[^f~l^+  .  .  . 

Only  the  first  two  terms  in  this  series  will  be  used  bee 
is  assumed  to  be  small.  'I'he  angle  0  is  to  be  determin 
the  minimum  mean-square  deviation  of  the  stream  fu 
from  its  mean  value.  Because  the  object  of  the  rotai 
essentially  to  adjust  the  shape  of  the  nose,  the  rotation 
also  bo  made  to  reduce  the  root-mean-squaro  deviation 
stream  function  to  a  minimum  for  a  portion  of  the 
including  the  nose. 

The  mean  value  of  the  stream  function  at  any  angi 

?(«-!/;  #(,.«  *-ij;  * 

The  diflercnce  between  the  new  stream  function  p  (s,  0 
its  mean  value  p  (0)  is  squared  and  integrated  to  obi 
measure  of  the  variation  of  p  (s,  0)  from  the  mean  vn 
the  now  angle.  The.  condition  for  obtaining  a  rniti 
root-mean-square  deviation  by  adjusting  0  is 


-I  {«».«) + 1>  -  ♦(«]  pr  -  ■“ 


The  second  integral  vanishes  by  virtue  of  equation  (11), 
which  may  also  be  used  to  eliminate  p  (0)  from  the  remaining 
term.  The  solution  for  0  is 


fVw* 
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In  order  to  apply  equation  (14),  dp/d0  must  be  known  at 
points  along  the  boundary  of  the  airfoil.  For  the  isolated 
airfoil,  the  contribution  of  the  vortices  is  unaffected  by  the 
rotation  and  therefore 


<lp_<H'u_  d 
d0~^  d0  ~il0 


(.-x\\+yV,) 


V  T/ 


(15) 


If  the  airfoil  is  rotated  about  the,  point  (fc,  i/c),  equation  (15) 
becomes 


^-|=cos  0  [(j— JV)  (/+(;/ -y/EJ-h 

sin  0  [{x-xc)  V,  —  {y—yP)  h/] 

where  (r,  y)  arc  the  coordinates  of  the  point  before  rot: 
For  small  values  of  0,  equation  (16)  reduces  to 

f-={x-x,)V.+  {y~i,,)V, 

The  choice  of  (x^,  y  j  will  have  no  effect  on  the  results  ii 
case. 

When  the  airfoil  in  cascade  is  rotated,  the  change  i 
position  of  the  vortices  of  the  adjacent  blade  must  be 
sidered.  For  the  isolated  airfoil,  it  was  unnecessary  ((> 
sider  the  change  in  position  of  the  vortices  becausi 
influence  of  a  vortex  (equations  (3),  (4),  and  (5))  depi 
on  the  fuuction/i,  which  is  constant  on  circles.  The  infli 
of  the  vortices  on  the  airfoil  is  therefore  independei 
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lliipcl ion.  Bi'ciuisc  ihc  /j  contours  iii'c  not  circles,  the  roiii- 
tion  in  cascade'  does  have  an  ed'oet,  which  is  approxiiuateil 
l)y  eonsiderin<j  all  closed  /2  curves  (/o<^0)  as  circles  in  order 
that  the  ell'c'ct  of  all  vortices  in  the  ri'^ion  /'..<()  may  he 
ne>rlecti'd  during  rotation.  The  ell'eet  of  all  vortices  in  the 
rcfrion  is  estimated  hy  assutniiii;  that  all  ihe/o  contours 
for  ./T-'O  are  straifilit  lines  uniformly  spaced.  The  How 
corres[K)nds  to  that  hi'tween  two  infinite  slraiirht  parallel 
vorte.x  sheets  of  uniform  stren<rth  tier  unit  leiifrlh.  This 
How  induced  by  the  vortices  in  the  re<rion  /2TO  is  in  the 
/-direction,  and  the  direction  of  the  How  indiiei'd  hy  the 
vortices  for  which  is  opposite  in  sense  10  that  induced 

by  the  vortices  for  which  .Vo<,i/. 

As  tlu'  point  hi'inij  considered  is  chan}/e<l.  the  rcfiions  for 
and  i/o<.V  will  include  diH'erenI  sections 

of  th('  blades,  and  hence  different  vortieity,  with  the  result 
that  the  /-velocity  component  c,  induced  hv  the  vorlex 
sheets  will  vary  with  the  point  under  consideration.  The 
alfrehraic  sum  of  the  /-component  of  the  imifonu  flow 
velocity  and  the  variable  /-velocity  fi  induced  hy  the  vortices 
in  the  region /.>()  is  to  he  used  like  the  velocity  component 
f’i  in  rotation  of  the  isolated  airfoil  (ocpiation  (17)).  The 
(luantity  ly  in  eipiation  (17)  is  replaced  by  the  corresponding 
I’j  ,=  I'l-fi',,  The  vortex  strength  per  unit  length  at  any 
point  on  the  airfoil  is  equal  to  ^'(So)  and,  therefore,  from 
equation  (10)  the  /-component  of  the  velocity  induced  by  the 

vortices  is  (p'{na)  duo,  where  the  integration  is  carried 

out  over  the  portion  of  the  airfoil  where  f>(s,$o)y-0.  A 
distinction  must  be  made  between  the  two  regions  .i/o<(,'/  'totl 
y„>i/  because  the  induced  velocity  components  liave  oppo¬ 
site  directions. 

The  computed  result  of  rotating  an  airfoil  in  cascade  de¬ 
pends  uiion  the  choice  of  (/eiVc)-  In  order  to  minimize  the 
error  involved,  values  of  d'pld&  are  reduced  by  choosing 
(/t,)/r)  as  the  centroid  of  the  vortex  distribution  on  the  airfoil. 
If  the  improvement  in  the  mean-square  deviation  of  ^  is 
small  compared  with  its  original  value,  it  may  he  preferable 
to  omit  the  rotation  of  the  airfoil  because  of  the  error  inherent 
in  the  approximation  for  dip jd0.  The  decision  should  he  made 
chiefly  on  how  ip  varies  at  the  airfoil  nose  and  whether  it  is 
approaching  a  constant  value  in  this  region  with  successive 
corrections  of  the  shape. 

Distortion  of  the  shape, — The  stream  function  computed 
after  the  isolated  airfoil  has  been  rotated  will,  in  general, 
still  vary  along  the  boundary.  This  variation  can  he  re¬ 
duced  hy  distorting  the  shape  of  the  airfoil.  If  the  distor- 
lion  is  small,  the  change  in  distance  between  any  two  points 
on  the  boundary  will  be  small,  although  the  change  in  the 
direction  of  a  segment  joining  those  points  may  he  consider- 
abh'.  The  eli'ect  of  the  distortion  on  the  contribution  to  . 
the  stream  function  of  the  vortices  on  the  boundary  is 
consequently  neglected.  The  largest  effect  of  the  distortion 
will  be  to  change  the  position  of  the  boundary  points  in  the 
uniform  stream.  The  airfoil  is  therefore  distorted  in  such  a 
manner  that  the  change  in  the  contribution  of  the  uniform 
•Stream  to  the  stream  function  will  eliminate  the  variations 
in  stream  function.  For  points  directly  opposite  each  other 
on  the  airfoil,  the  change  in  distance  will  be  of  the  same  order 
of  miignitudc  as  the  distortion.  Consequently,  distortions 


that  result  in  change  of  thickness  of  the  airfoil  converge 
very  slowly  heeause  of  the  inaiteuracy  of  the  ftnulamental 
assumption  on  which  the  correction  is  hasetl. 

'I’htis,  when  the  stream  ftmetion  along  the  hotmdary  of 
the  isolated  airfoil  is  known,  some  luimbcr  is  arbitrarily 
ehoseti  as  the  desired  cotistant  value,  of  the  stream  function. 
If  A\p=ip  —  ip  is  the  tliH'eretice  between  the  comptited  stream 
ftmetion  tit  a  jioint  atul  the  desired  cotistant,  the  point  is 
moved  a  distance  —  Ai/'/l'  perpendicular  to  the  direction 
of  the  mean  stream,  where  the  direction  of  increasing  uni¬ 
form  stream  function  is  taken  as  positive.  The  airfoil  in  a 
cascade  is  distorted  in  tlu'  same  maimer,  by  using  the  varying 
resultant  local  mean  stream  velocity  -\/T’i.r'-f- 1 V;  corrections 
are  made  with  ip  equal  to  the  mean  value  of  \p  on  the  airfoil. 

COMPUTATIONAL  PROCEDURE  FOR  CASCADES 

CHOICE  OF  VELOCITY  DISTRIBUTION 

.'several  factors  iidluenee  the  choice  of  the  velocity  distri¬ 
bution  for  which  an  airfoil  is  to  be  found.  Especially  in 
rotors,  sturdy  blades  are  required.  Long  thin  tail  sections 
must  be  avoided  and  where  high  rotative  speeds  and  stresses 
occur,  overhang  of  thin  sections  is  likely  to  induce  blade 
failure.  The  radial  distribution  of  airfoil  cross-sectional 
area  is  also  fundamental  in  determining  the  blade-root 
stresses.  Overhang  can  be  reduced  by  proper  choice  of  the 
velocity  diagrams  for  the  sections,  but  the  other  factors  are 
influenced  chiefly  by  the  thickness  of  the  section. 

'I'he  desired  thickness  may  be  attained  by  first  assuming  a 
blade  shape  and  spacing  and  by  then  using  the  stream-filament 
method  of  reference  6  to  compute  the  velocity  distribution 
over  a  portion  of  the  airfoil  that  determines  the  thickness. 
The  spacing  may  be  regarded  as  fi.xed  but  the  curvature  can 
be  adjusted  if  local  velocities  are  too  high  for  the  desired 
thickness.  This  computed  velocity  will  then  serve  as  a 
guide  to  the  choice  of  an  airfoil  velocity  distribution,  which 
should  be  chosen  to  avoid  high  velocity  peaks  and  steep 
negative  gradients.  If  the  average  of  the  velocities  on 
opposite  sides  of  the  blade  camber  line  is  retained  in  the 
modilieation  of  the  velocity  distribution  computed  from  the 
stream-filament  method,  the  thickness  will  also  be  retained. 

Because  of  the  irrotationality  of  the  fluid  motion,  the 
velocity  integral  or  circulation  around  the  airfoil  must  be 
eijual  to  that  around  a  blade  but  over  a  width  equal  to  one 
blade  space.  Therefore, 

,  J  ¥>'(■''■)  d'>'=r=>S’(F2.,i  — T-Vd 

where 

r  circulation  about  airfoil 

I'r.i  tangential  velocity  entering  cascade 

l'j.,2  tangential  velocity  leaving  cascade 

This  relation  places  a  condition  on  the  assumed  velocity 

distribution. 

If  the  computations  thus  far  have  been  made  in  order  to 
select  a  velocity  distribution  for  the  airfoil  cascade  in  a  com¬ 
pressible  fluid  flow,  an  eipiivalent  velocity  distribution  for 
the  flow  of  an  incompressible  fluid  must  be  determined 
before  the  blade  shape  can  be  computed  by  any  method 
based  on  incompressible-flow'  theory.  For  stibcritical  flows, 
the  directions  of  the  inflow  and  discharge  velocities  are 
nearly  the  same  for  compressible  and  incompressible  flows, 


f) 


HKPOUT  XO.  SCiil - XATIOXAL  AD\  IS(>i;V  COMMITTEE  FOR  AEHOX AF'IICS 


bill  for  iiu'oinpi'('ssil)lo  How  the  coinpoiii'iil  nonnnl  to  casciKlc 
axis  is  the  same  upslream  aiici  downslri'ani.  The  Ivarmiiii- 
Tsien  compressibility  corn'clioii  ( rcforciicc  7)  or  tliat  of 
Garrick  and  Kaplan  (reference  Sf  may  be  ai)])lied  to  the 
velocity  on  tbe  blade  surface  to  estimate  roughly  the  corres-  j 
pondiny  iucompressibb'-llow  velocity  distribution.  Tbe  re¬ 
sultin';  V('locity  distribution  in  any  ease  must  satisfy  the  i 
circulation  condition.  This  |)rocedure  does  not  t;ive  an  | 
exact  solution  lor  com])ressible  Hows,  but  the  resultant  com- 
pressibb'  How  will  have  a])proximately  I  be  desired  ebar- 
acferistics  of  low  ])r('ssure  Gradients  and  no  bisrb  velocity 
peaks. 

COMPUTATION  OP  AIRFOIL  SHAPE  FROM  THE  CHOSEN  VEI.OCITV 
OLSTRIIIUTION 

The  numerical  computation  of  tbe  (|uanlitie.s  involved  in 
tbe.  preceding  analysis,  particularly  tbe  function  /-,  is  ex-  j 
tremely  laborious  wben  tables  ol  are  used.  Most  of 

tbe  computations  are  tborefore  executed  grapbically.  In  tin- 
cascade  example,  tbo  air  was  assumed  to  enter  tbe  cascade 
at  an  angle  of  45°  from  tbe  cascade  axis  and  to  leave  at  an 
angle  of  —110°  from  tbe  cascade  axis.  Tbe  prescribed  veloc¬ 
ity  distribution  is  given  in  figure  Hfa).  Tbe  value  of  tbe  lift 
coefficient  for  this  airfoil  is  3.1.  Tbe  sbapes  of  tbe  isolated 
airfoil  and  tbe  airfoil  in  cascade  are  comptited  by  tbe  fol¬ 
lowing  stejis: 


(:0  Initial  nirfnil. 
(h)  Final  airfoil. 


FiGUnE  :i.— Trcscriiicti  rolocitv  tiislribiilion  for  Ihiok  airfoil  In  casoiuli*.  t 

j 

1.  Curves  for  constant  /,  for  tbe  isolated  airfoil,  or  con-  ! 
stant/a  (lig.  2)  for  tbe  airfoil  in  cascade,  are  drawn.  Tins  j 
diagram  sbould  be  made  on  some  transparent  material  that  ; 
will  ebange  neither  in  size  nor  shape.  Tbe  coordinates  of  tbe  i 
curves  for  constant /j  are  given  in  table  I. 

2.  A  desired  velocity  ip'{g)  is  chosen  as  a  function  of  tbe  i 

arc  length  of  tbe  airfoil  (fig.  ;i(a)).  An  airfoil  sba])c  having  i 
the  desired  total  arc  length  is  assumed  and  is  draAvn  to  tbo 
same  scale  as  the  plot  of/,  ov  f 2-  Tbe  draw'ing  is  made  on 
grid  paper  and,  in  the  case  of  a  cascade,  the  ai-axis  coincides 
with  the  cascade  axis  (fig.  4).  ! 

3.  The  velocity  distribution  ip'(s)  is  integrated  to  obtain  I 
(he  velocity  potential  <,r(.'-f.  This  function  is  plotted  on  tbe  ! 


same  chart  as  tin-  assumcil  airfoil  sbapi-  foi'  the  cori'cspi 
//-coordinate,  as  shown  in  figure  4.  by  plotting  both 


rioVRE  4.— I’lot  of  airfoil  and  volocit.v  potential  for  use  in  coinputatitm. 


the  //-coordinate  of  tbe  airfoil  against  a-  on  a  supplcnn 
graph.  In  regions  of  the  airfoil  where  //  varies  little  v 
(bat  is,  where  tbcairfoil  boundary  is  parallel  totbcj-dirc 
a  sbould  be  plotted  against  x  in  tbe  same  manner,  as  s 
in  figure  4. 

4.  In  ord(>r  to  find  tbe  stream  ftmetion  at  a  point 
on  tbe  airfoil./)  must  be  plotted  as  a  ftmetion  of 
to  evaluate  (be  rpiantity  J'f2(K,Sa)  dip(s„)  of  ccpiatioi 
If  tbe  chart  of/j  is  suiterimposed  on  tbo  airfoil  with  om 
tex  center  overlying  tbe  ])oint  (x,?/),  tbe  value  of  /^  m. 
rend  at  (x„,,i/„)  and  tbe  corresponding  value  of 
also  1)1-  read  from  tbe  plot  of  ip{Xo,yc).  Tbe  value  of  / 
is  tbe  same  as  would  have  been  obtained  by  centerin 
chart  on  (x„,,i/„)  because  of  tbe  symmetry  of  tbe  fun 
A  succession  of  values  of  is>  and/s  are  obtained  in  this  fa 
for  various  positions  (j’„,//„)  that  intersect  tbe  ,/>  con' 
and  a  ])iot  of  these  iioints  (/2,'p)  may  be  made  for  tbe  ass 
|)Osition  (x.y).  This  /troccdure  is  illustrated  in  figure 
a  itarticular  point  (x.y)  on  which  tbo  f>  chart  is  ccni 
Tbe  readings  for  a  jtarticular  {x„.yo)  are  shown  by  the  an 
lines.  The  points  L  to  (i  on  tbe  blade  are  shown  01 
corresponding  /j  curve.  Tbo  discontinuity  of  ip  bet 
points  1  and  0  is  tbe  circulation.  Tbe  iliscontinuity  be: 
4  and  5  represents  tbe  region  where  /.  approaches  —  o' 

.5.  The  proper  method  of  integration  then  proceeds 
I  through  ()  to  7  and  then  to  (be  origin,  with  consti 
front  4  to  5.  Tbe  region  from  4  to  5  with  /j  appi'oa 
—  00  is  computed  by  e(|uation  (ti)  or  (Ob);  the  coustan 
assumed  to  be  tbe  radius  of  the  near-circle,  wbicb  c 
ponds  to  tbe  value  of  /  where  tbe  discontinuity  from 
.5  occurs. 

The  total  area  including  this  small  addition  is 

J  ip/Sol/As.So)  (/s= J'/2(a’,.s„)  dip 

wbicb  is  tbe  stream  function  due  to  vortices  on  the  entii 
of  airfoils  in  cascade.  IMiere  /2=0  at  tbe  points  .1,  1 
and  D  (fig.  5),  tbe  values  of  p  arc  noted  as  ».i(.s),  pb(-'^).  < 
and  pnid).  These  values  arc  used  in  computing  tbe  str 
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function  change  caused  by  rotating  the  blade.  The  streani 
function  at  the  jioint  {x,  ij)  may  now  be  computed  from 
equation  (8)  or  (o),  and 

—  f  >-f 


Fiqure  ^.—Superposition  o(  Hgurcs  2aod  4  to  obtain  plot  of /s  nsainst 

A  plot  of  the  stream  function  (variation  from  the  mean 
value)  is  shown  in  figure  G  for  the  initially  assumed  shape. 
Corresponding  points  on  adjacent  airfoils  have  a  difference 
of  A\l//ViiS  equal  to  1.0. 


0  Initial  alrfo/l\  \  I 
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Figure  0.— Variation  In  stream  function  alouK  initial  shape  and  first  and  seventh  appro.xi- 
matlons  of  airfoil  cascade. 

G.  Wlien  \p(s)  is  known  at  a  sufficient  number  of  points, 
the  airfoil  may  be  rotated  as  previously  described.  For  the 
isolated  airfoil,  equations  (14)  and  (17)  may  be  used  directly. 
For  the  airfoil  in  cascade,  the  coordinates  of  the  centroid  of 
the  airfoil  must  first  be  computed  by 


< 

^xp'is„)  ds„ 

Before  equation  (17)  can  bo  useil  to  compute  the  vari¬ 

able  quantity  F,,,  must  be  computed.  The  vortices  in  the 
region arc  considered  to  be  uniformly  distributed  along 
the  cascade  axis  and  the  velocity  induced  by  such  a  distri¬ 
bution  is 

±  % 

where  7  is  the  vortex  strength  per  unit  length  along  the  cas¬ 
cade  axis  for/2)>0.  Therefore, 


where  tlio  integral  is  to  be  taken  o\Tr  the  regions /2>0.  The 
region /2>0,  ?/o>i/  contributed  a  positive  component  to  v^, 
whereas  the  region  //„<;/  contributes  a  negative  com¬ 

ponent.  The  computation  is  simply  carried  out  by  making 
use  of  the  fact  that  the  integral  for  Vx  is  the  difference  between 
values  of  (p  at  points  where  fi—O.  The  values  of  <pa.{So), 
‘PaiSo),  yc(«o),  and  v’afSo)  from  step  5  are  used  at  this  point 
to  obtain 

2t>xS=  yifi'lso)  dSo~<PA  —  yo-[-r — (<pc — ipa)  (18) 

where  T  is  introduced  because  of  the  discontinuity  in  p  at 
the  trailing  edge.  The  sum  (Pa—vd+T  gives  the  effect  of 
the  vorticity  in  the  region,/2(s,  s„)>0  near  the  trailing  edge, 
and  the  term  0c— pb  gives  the  effect  of  the  vorticity  in  the 
region/2(s,s„)>0  near  tlio  leading  edge.  If  either  the  leading 
edge  or  the  trailing  edge  lies  in  the  region /2(«,s„)<0,  only 
two  points  of  intersection  will  remain  and  one  of  the  two 
groups  of  terms  in  equation  (18)  will  vanish.  The  quantify 

—  j*  p'{s,)  ds„  is  added  to  the  x-component  of  the  original 

uniform  stream  velocity  and  the  quantity  d>l//d0  of  equa¬ 
tion  (17)  may  be  computed  for  a  number  of  points  and  the 
angle  0  computed  from  equation  (14),  using  the  values  of  {XcPjc) 
just  determined.  After  these  computations  have  been  made, 
the  airfoil  is  rotated  through  the  angle  0,  and  the  value 

. 

IS  assigned  as  the  value  of  the  stream  function  of 
the  point  after  rotation. 

7.  A  value  of  is  known  at  points- along  the  airfoil 
boundary.  The  mean  value  over  the  airfoil  ^  is  subtracted 
from  ^  loaving,A^.  For  the  isolated  airfoil,  no  subtraction 

is  necessary.  Each  point  is  moved  a  distance - 7== 

^|VxJ+l7 

in  the  direction  perpendicular  to  the  velocity  computed  in 
step  6.  The  curve  joining  the  points  in  their  new  positions 
is  the  adjusted  airfoil. 

8.  The  total  arc  length  of  the  adjusted  airfoil  will  be 
different  from  the  original  one,  in  general,  although  local 
changes  in  length  ivill  be  negligible.  The  airfoil  is  so  scaled 
that  the  length  of  the  suction  side  is  the  same  length  as  it  was 
before  distortion  because  this  surface  is  the  critical  surface 
of  the  airfoil.  This  process  will  result  in  a  change  in  length 
of  the  pressure  side.  The  velocity  over  the  pressure  side  p'  (s) 


8 


REPORT  NO.  SG9 — X.\TIO.\AL  .ADVISORY  COMMITTEE  FOR  AEROX.AUTICS 


imist  then  he  altered  in  such  a  manner  tliat  tlie  difrerencc 
in  potential  between  the  two  stagnation  points  remains  the 
same.  .Vs  a  result,  the  quantities  that  retain  specified  values 
are  the  lensrth  and  the  velocity  distribution  on  the  suction 
side  and  the  circulation  around  the  airfoil.  The  entire 
procedure  is  repeated  with  the  adjusted  shape  until  the 
variations  in  the  stream  function  result  in  very  little  chansrc 
in  the  shaiie  of  the  airfoil. 

mSCU.SSION  OF  EXAMFI.ES  AND  TECHNIOUES 

For  the  e.vamplc  being  computed,  the  stream  functions 
obtained  for  the  initially  assumed  shape  and  the  first  and 
seventh  aiipro.ximalions  are  plotted  against  the  arc,  length 
(fig.  0),  which  is  taken  as  zero  at  the  trailing  edge  and  pro¬ 
ceeds  counterclockwise  around  the  airfoil  as  shown  in 
figure  7.  The  fact  that  Atp  for  the  initial  shape  is  ])ositive 
over  the  first  half  of  the  arc  length  and  negative  over  the 
second  half  indicates  that  it  is  too  thick  because  the  required 
distortion  in  shape  will  make  it  thinner.  The  change  in 
thickness  results  in  a  change  in  velocity  distribution  over  the 
pressure  side  of  the  airfoil  in  order  to  maintain  the  desired 
circulation.  The  velocity  that  was  originally  assumed, 
which  is  equal  to  the  vorticity  per  unit  length  distribtited  on 

......  Inifial  Qirfoil 

Pinal  airfoil 


Figure  7.— Initial  shspo  and  Anal  approximation  of  thick  airfoil  showioR  cascade  spacinR. 

the  initial  airfoil,  is  shown  in  figure  3  (a)  and  the  velocity 
over  the  final  shape  in  figure  3  (b).  The  length  of  the 
pressure  side  has  increased  and  the  velocity  has  decreased 
in  the  portion  of  1:1.1. 


Over  the  section  of  the,  airfoil  that  has  collapsed  ' 
thickness,  the  surface  velocities  of  figure  3  (b)  may  m 
been  obtained,  but  the  loading  (circulation  per  m 
length),  which  is  the  difference  in  the  velocities  on  oi 
sides,  has  been  realized.  In  practice,  this  collapse 
vented  by  increasing  the  assumed  velocity  on  the 
surface. 

If  the  initially  assumed  airfoil  shape  has  a  thickne 
differs  considerably  from  the  correct  one,  the  process  oi 
adjustment  will  converge  rather  slowly.  The  labor 
reduced,  however,  by  computing  the  stream  fiinctio 
few  points  on  the  airfoil  and  locating  these  points  to 
mine  the  thickness.  This  procedure  is  followed  for  tl 
few  appro.xirnations  until  the  thickness  of  the  airfoil  i^ 
accurate.  The  stream  function  is  then  computed 
larger  number  of  points,  particularly  near  the  leadin' 
in  order  to  get  more  detail  of  the  shape. 

Arbitrary  specification  of  a  velocity  distributioi 
result,  not  in  a  physically  real  airfoil,  but  in  a  1 
shape  or  a  collapsed  shape  (zero  thickness  over  a  por 
the  blade).  The  velocity  distribution  must  then  be  m 
to  obtain  a  real  shape;  these  modifications  should  be  si 
to  keep  the  desirable  properties  of  the  original  distril 
Velocity  peiiks  and  steep  velocity  gradients,  which  t, 
occur  on  the  suction  side  of  an  airfoil,  are  to  be  ax 
If  the  airfoil  collapses,  the  vorticities  of  the  two  sides  t 
cancel  each  other  and  the  remaining  vorticity  rcpresei 
difference  in  velocity  across  the  thin  airfoil  rather  th, 
velocity  along  the  boundary. 

The  method  was  also  applied  to  the  design  of  a  thin 
(camber  line)  in  a  cascade.  The  vortex  distribui 
equivalent  to  load  distribution  (difference  in  velocity 
the  airfoil)  rather  than  velocity  as  in  the  case  of  a 
airfoil.  The  velocity  diagram  for  the  cascade  and  the  ( 
load  distribution  for  the  thin  airfoil  are  shown  in  fig 
The  value  of  the  lift  coefficient  of  the  resultant  airfoil 


Figure  8.— Velocity  diagram  for  cascade  and  prescribed  load  distribution  for  thin 

cascade. 


ISOLATED  AKD  CASCADE  AIRFOILS  WITH  PRESCRIBED  VELOCITY  DISTRIBUTIOX 


9 


The  iniiiiil  shape  was  ohiainod  by  assumine:  zero  spacinc; 
between  the  airfoils.  The  initial  shape  and  the  (irst  and 
third  a[)pro.\iinations  to  the  airfoil  shape  are  shown  in  lisrure  il. 

. Initial  oir-f dll 

- First  approximation 

-  Third  approximation 


< - s 

Figure  9.— Assumed  shape  aod  first  and  third  approximations  of  thin  airfoil  showing  cascade 
spacing. 

The  second  and  third  approximations  differ  very  little. 
The  third  approximation  is  redrawn  in  tliis  diagram  to  show 


the  spacing  between  airfoils.  The  convergence  of  the  method 
is  shown  graphically  in  figure  10.  The  variation  Ai/-  of  the 


Figure  lO. — Variation  in  stream  fiiiiclion  for  successive  approximation  of  thin  airfoil  in 

cascade. 


stream  function  from  its  mean  is  divided  by  VS  to  make  it 
dimensionless  and  is  plotted  against  the  arc  length  along  the 
airfoil  whore  s=0  at  the  trailing  edge.  The  stream  function 
computed  on  the  second  approximation  is  nearly  constant, 
which  gives  the  thirrl  approximation  almost  the  same  shape 
as  the  second  one.  The  rapid  adjustment  of  camber  con¬ 
trasts  with  the  slow  adjustment  of  thickness. 


Flight  Propulsion  Research  Laboratory, 
National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  March  4,  19^7. 


APPENDIX  A 

VELOCITY  DISTRIBUTION  ON  THE  DERIVED  AIRFOIL  AT  DIFFERENT  FLOW  ANGLES 


Conformal  mapping. — Wlieu  an  airfoil  is  given,  the  ve¬ 
locity  distribution  over  its  surface  must  frequently  be  found 
at  different  angles  of  attack.  This  problem  may  be  solved 
by  tbo  method  of  conformal  mapping,  which  consists  in 
mapping  the  region  exterior  to  the  airfoil  on  the  exterior  of 
a  circle.  The  velocity  around  the  airfoil  is  obtained  from  the 
known  velocity  around  the  circle.  Procedures  for  finding 
the  function  that  maps  a  given  airfoil  into  a  circle  are  pre¬ 
sented  in  references  1  and  9  for  the  isolated  airfoil  and  refer¬ 
ences  ;i  and  10  for  the  airfoil  in  cascade. 

In  general,  the  procedure  for  finding  the  mapping  function 
of  an  airfoil  is  a  laborious  one.  But  when,  as  in  the  present 
case,  the  velocity  distribution  over  the  airfoil  at  a  particular 
angle  of  attack  is  known,  the  correspondence  between  points 
on  the  airfoil  and  on  the  circle,  and  hence  the  How  velocity 
at  other  angles  of  attack,  can  be  obtained  very  easily. 
Indeed,  the  correspondence  of  points  and  the  velocities  for 
various  angles  of  attack  can  bo  obtained  by  the  method  given 
in  reference  11  from  the  initial  data  without  knowing  the 
airfoil  shape,  because  the  complex  potentials  of  the  airfoil 
plane  and  the  mapping-circle  plane  are  equal.  Before  the 
airfoil  is  designed  it  is  therefore  possible  to  check  whether 
the  airfoil  to  be  computed  will  be  satisfaetory  under  condi¬ 
tions  diflercnt  from  the  design  condition. 

Isolated  airfoil. — The  How  about  any  airfoil  shape  can  be 
mapped  on  the  How  about  a  unit  circle  in  such  a  way  that 
corresponding  points  have  the  same  potential.  The  How 
about  the  airfoil  is  given  and  the  potential  function  at 
each  point  is  computed.  If  the  potential  function  on  the 
airfoil  is  computed  by  integrating  the  velocitj'  from  the 
stagnation  point  at  the  trailing  edge  in  a  counterclockwise 
direction  around  the  airfoil  oriented  as  in  figure  1.  the  poten¬ 
tial  will  be  zero  at  the  trailing  e<lge.  decrease  to  a  minimum 
<Pmin  at  the  stagnation  point  at  the  leading  edge,  and  then 
increase  to  a  value  equal  to  the  circulation  T  at  the  trailing 
edge.  The  corresponding  flow  about  the  circle  is  determined 
by  the  conditions  that  it  must  have  the  same  values  of 
and  r  for  a  cotTcspondence  to  exist  between  all  airfoil  and 
circle  points.  If  8t  is  the  central  angle  of  the  stagnation 
point  on  the  circle  that  corresponds  to  the  trailing  edge  of 
the  airfoil, 

^!^=-(cot0r+9r+W2)  (Al) 


Equation  (Al)  can  be  solved  numerically  for  8t  because  all 
the  other  quantities  are  known.  The  velocity  at  infinity  in 
the  circle  plane  Vc  can  then  be  dcteiTnined  from  the  Kutta- 
Joukowsky  condition,  which  requires  that  Sr  be  a  stagnation 
point;  that  is, 


The  velocity  potential  at  points  on  the  circle  is 

pa  p 

(Pe~—2Vc  COS  0+o''  cos 

-TT  Jtt 

p 

The  quantity  2  Vc  cos  dr—  Or  is  a  constant  that  is  suhtr 

-TT 

ill  order  to  make  at  the  .stagnation  point  corresi 

ing  to  the  trailing  edge. 

The  correspondence  of  points  on  the  airfoil  with  poiu 
the  circle  is  obtained  by  associating  points  where 
The  velocity  on  the  circle  at  a  uniform  stream  flow  angl 

Vc(8,a)  =  2Vc  [sin  (9-fa)— sin  (Sr-f  a)] 

The  nature  of  the  conformal  transformation  is  such  thii 
ratio  of  the  velocity  at  a  point  on  the  airfoil  to  the  vel 
at  the  corresponding  point  on  the  circle  is  independc. 
angle  of  attack.  Therefore,  the  velocity  «?„'(■'*)  on  the  ii 
at  flow  angle  a  is 

Va'i^)  __  y'(*’) 

Vc{e,a)  L'c(e,o) 

where  the  design  flow  angle  is  taken  as  zero.  Equation 
can  be  used  to  compute  the  velocity  distribution  on  tin 
foil  e.xcept  at  the  two  points  that  were  stagnation  poin 
the  design  angle  of  attack. 

Airfoils  in  cascade. — The  flow  about  a  cascade  of  aii 
can  be  mapped  conformally  into  the  flow  about  a  unit  i 
with  two  singular  points  located  on  the  real  axis  symn 
cally  with  respect  to  the  center  of  the  circle.  These  sin; 
points  correspond  to  the  points  at  infinity  in  front  of 
behind  the  cascade,  respectively.  In  a  cascade  of  air 
the  distance  of  these  points  from  the  center  of  the  cin 
uniquely  determined  by  the  same  conditions  that  dctci- 
the  flow  about  the  circle  in  the  isolated  case;  namely 
circulation  per  airfoil,  the  velocity  potential  at  the  leu 
edge,  the  blade  spacing,  and  the  upstream  and  downsti 
flow  angles. 

The  distance  from  the  singular  points  to  the  center  o 
circle  is  denoted  by  e‘^.  The  flow  about  the  circle  is 
that  the  location  of  the  stagnation  points  8,  is  doterm 
by  the  relation 

r  sin  8,  ^  ,  cos  8,  .  , 

2ES““sinh  ^  ^ 

where  X  is  the  angle  of  inclination  of  the  mean  stream  t( 
normal  to  tbo  cascade  axis.  (See  reference  0  for  dcti 
The  quantities  T,  V,  S,  and  X  are  known  from  the  flo' 
the  cascade  plane  and  therefore  equation  (A6)  provid 
relation  between  K  and  the  location  of  the  stagnation  pm 
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'I’iu'  vcldcily  potcMtiiil  iit  any  poiiil  on  llio  circle  is 

US’/  _  sin  0  ^  ,  1  -1  ^  ^  I  I’  ,  -1  ^ 

(  Sin  X  liin  '  ■  ,  r-— cos  X  tiiiih  ‘  t— - tan  ‘  , 

TT  \  sinh  K  cosh  A/  tanli  A 

r'  *''/  ■  ^  ,  -1  sin  0,-  ,  I  cos  07- \  I  r  _  tan07.“| 

X  tan  ^^--cos  X  tank  ,,,,,1^7,-)  +  .,^  tan  (Ai) 

where  O,-  is  the  particular  value  ot  0^  corresponding:  to  the  trailinir  ed>rc  of  tlu'  airfoil.  Tlie  e.vpression  in  hi’aekets  is  a 
eonsiant  so  chosen  that  the  |)otential  will  vanish  at  the  stagnation  point  coi-responding  to  the  trailing;  edge  of  the  airfoil.  In 
orih'r  to  map  the  cascade  on  the  circle,  K  must  he  found  so  that  the  value  of  at  fly,  the  stagnation  point  6,  eorresjionding 
to  till'  leading  I'dgi'  of  the  airfoil,  is  e(|nal  to  the  value  of  the  velocity  [lotential  there.  The  identity 


(sS  a)'  ^+(cS4:)’ 


is  used  to  eliminate  j-.from  ecpiation  (Ati)  to  give 

sin  fl,  cosh- A' sin  X±cos  X-^cosh*  A— cos- X— cosh*  A  sinlr  A 

sinh  A  coslr  A— cos-  X 

In  successive  appro.ximations,  a  value  of  K  is  assumed  and  equations  (AS)  and  (AG)  are  used  to  find 

Jfiili  A'’  c*osfi^  These  values  are  inserted  into  equation  (A7)  to  find  ipc.e  al  6=0.v.  If  yc.c(9.v)  is  not  equal  to  ipmin- 

another  value  of  A  is  chosen,  on  the  premise  that  <Pc.ci?it)  "‘ill  decrease  as  A  is  decreased.  When  is  evaluated, 

care  should  he  taken  to  use  consistent  values  of  the  inverse  tangents.  After  two  values  of  A  and  (PeA^.v)  are  determined, 
interpolation  or  e.xtrapohition  may  he  used  for  new  values  ot  A, 

When  A  has  been  found,  it  is  used  in  equation  (A7)  to  evaluate  pe.c  at  values  of  fl  all  around  the  circle,  A  point  on 
the  circle  corresponds  to  the  point  on  the  airfoil  where  <p(i/)=ipex-  The  velocity  at  the  point  fl  on  the  circle  is 


_VS  sinh  2  A _  r  /  cos  fl  cos  flr  \  .  /si 

cosh  2A— cos  2fl  L  ^  \cosh  A  cosh  A/’^*’''^  Vsin 


■_sin_fl 
sinh  A 


sin  0T  ' 
sinh  Rj 


and  the  v-elocity  <,5a^(s')  on  the  airfoil  at  any  other  mean 
flow  angle  X  +  a  is 

V>a'Os')  =  'V.c(fl,X  +  ai  -;'g-xj  (AlO) 

as  in  the  case  of  the  isolated  airfoil. 

The  designed  airfoil  was  mapped  on  the  unit  circle  by  the 
method  described.  The  constant  A',  the  natural  logarithm 
of  the  distance  from  the  singular  points  to  the  center  of  the 
unit  circle,  is  0.075.  The  correspondence  between  points 
on  the  airfoil  and  those  on  the  circle  is  plotted  in  figure  II, 
which  shows  the  arc  length  of  the  airfoil  as  a  function  of  the 
central  angle  of  the  circle.  The  velocity  at  any  point  on  the 
airfoil  for  any  angle  of  attack  a  may  bo  obtained  from 
equations  (A9)  and  (AlO),  the  velocity  distribution  as  in 
(igure  3  (b),  and  the  relation  between  a  and  fl  as  in  figure  1 1. 

The  ratio  is  equal  to  (/fl/iA  (radians)  and  need  bo 

computed  only  once  for  any  given  airfoil. 


-200  -ISO  -mo  -50  0  50  100  150  200 

Central  angle  of  circle,  6.  degrees 

Figure  11.— Corrospoodenco  between  points  on  nirfoil  and  points  on  unit  circle  by  conformal 
transformation. 


APPENDIX  B 

DERIVATION  OF  THE  CASCADE  EQUATION 


An  equation  is  to  be  developed  for  tlie  complex  velocity 
at  any  point  in  the  Held  of  How  of  a  fluid  jtast  a  row  of  equally 
spaced,  confiruent  bodies.  Coordinates  a.xes arc  chosen  with 
the  origin  inside  one  of  the  bodies  and  the  x-axis  in  the 
tlirection  of  the  row.  (See  fig.  12.)  The  body  containing 


Figure  12.— Dlazram  for  derivation  of  equation  for  flow  about  ca.sca<lo. 


the  origin  is  denoted  by  /4,  bodies  along  the  posit 
rection  of  the  j:-axis  by  li\,  B2,  ■  ■  ■,  and  along  the  tu 
direction  of  the  r-axis  by  5_i, .  .  ..  A  circle  ^lo 
raditis  is  drawn  about  the  point  s  where  the  vclocit 
be  determined.  A  rectangle  R  is  drawn  with  its  cei 
the.  origin  and  its  sides  parallel  to  the  a.xes  of  length  (2.' 
and  width  2t,  which  contains  the  bodies  .  . 

Ro,  Rt-  ■  ■  ■  /?jv,  and  the  circle  A.  If  a  side  of  the  re( 
intersects  one  of  the  bodies,  the  side  may  bo  distorter 
around  the  body  with  no  e.ssential  change  in  the 
The  function  w'{z„)lz„—s)  is  an  analytical  function  ■ 
the  region  inside  the  rectangle  R  btit  outside  the  bod 
and  the  circle  A. 

Therefore 


JR^o  ^  Ja^^o  •* 


The  first  integral  can  be  broken  up  into  four  int' 
one  along  each  side  of  the  rectangle,  namely, 


Cw'iz,),  _  p 

Jb  2,-2  J  _(.V+,,2,S  + J-, 


w'\(N+\!2)S±iy, 


(N+\l2)S+iy-z 


Jw- 


t  .u-'ttor 
+i/2)s  x„-rtr— 4 


I 


w'[-(Af+ 1/2)5+ H. 
(A^+l/2)S+tt/,- 


In  an  evaluation  of  these  integrals,  the  function  w'{zA  is 
periodic,  with  period  5,  and  approaches  a  constant  value 
infinitely- far  from  the  cascade;  that  is, 


and 


'w'(x,-iij^)^W2'  as  Vo- 
re'(x„+/y„)-^w/  as  y„-^ 


From  the  last  of  these  conditions,  it  follows  that 
w'(x„— iV)=ie3'(x„— t'O  +  Wi' 

where 

tt’s'fxo— d)— >0  as  t— >00 

Therefore,  the  first  integral  on  the  right  side  of  eqi 
(B2)  is 


r(N+ims  w'(x„_,7)  ,  /•(.v+i.-2k'.'  ,/j.  j-(w-t-i,'2)s  '(j-  ;,) 

- =5 - (/x„=tt>,  - .: - h  - - “  «X„ 

./ -(A'+|.'2)S  X„ — ?f  —  2  J-(.V+1,'21SX„— tt  — 2  J_(.v+l/2)S  J„—lt—Z 


The  first  of  these  integrals  is 


^UV+I,'2)S  ,Ij.  [(,V-Li/2)S-d 


U-z] 


*inw^ 


as  and  03  ,  provided  that  t /{NS)  — ^t).  The  last  integral  in  equation  (BH)  is 

p.v+ 


■(.V+1/2)S  W3'(x„-d) 
.V+1/2)S  x„ — li — 2 


"■W 

II 

r("+''-'^'’'w3'(x„-d) 

dXg 

H  =  -.V 

J(n-I/2)S  X„ — it  —  2 

II 

>. 

w/{x„-it) 

dxo 

.J_.v/2  Xo+zyiS’ — it  —  ^ 

rsi2 

"^3  {Xg  ft)  J  , 

Xq - ll - H=1 

rsj2 

J-S12 

J  -  s/2 

A'  2(x„—it  —  z)w3'ix„—it)  dx„ 

(,c„—it  —  z)-—n~S'^ 
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If  t  is  cliospii  suiiiciciitly  largo  so  i  <e,  whore  t  is  any  prcassigned  positive  number,  the  moduli  of  the  integrals 

are  less  than 

rc.s'2  ,ix„  '2\x-it-z\dx„  n  r  r»-  dx„ _  ^  2  V(jT=^)-T'(f+^^  dxol 

‘  L)-s,2  x„—ii  —  z\^^i  J-K/2  \lxo—il—z)-—irS'-lJ  -  LJ-s'2  ^  (.j:,,— /)'+  (t+l/)-  n  =  i  Uo— •<■)■+  {t  +  ij)-  —  n-S-J 

When  A’-^oo  ,  this  quantity  approaches 

ej  cot  +  dx„ 


This  integral  is  finite  and.  because  «  can  be  made  arbitrarily 
small  as  0° ,  the  last  integral  in  equation  (B3)  approaches 
zero.  Therefore, 


.c: 


(•v+i/2)s  w'(r„-d) 


dXo-^Tritc,' 


as  l-^co  and  v  I'l  same  wav  and  under  the  same 

A'o 

conditions. 


j;; 


-(A-+1/21S  ir'(x„  +  !t) 


A- — ^  dx„-^inW2 
(A'+i;2)S  X„+lt—Z 

The  second  and  fourth  integrals  on  the  right  side  of  equa¬ 
tion  (B2)  can  be  evaluated  by  combining  them.  Because 
w'  is  periodic, 

w'  [  (A^+ 1  /2)  (■;/„]  =  w'  [  -  ( A^+ 1  /2)  5'+  f  y  J 

and  therefore, 


(,V+1/2)S  Xo—lt—S 

n  if'f(A^-bl/2)SJ-ii/„]  . ,  ,  f-' w'f-(A'  +  l/2)S+/i/„]  .J„  _  f  -2(A^-h  l/2)6Y’'[(A'-b  l/2)A^-b?.yj 

-(X+mS+ui7-z  -(N+ll2)S+iy.-z  J-,  (iy.-zy-(N+ 1/2^8^ 

The  velocity  w'[(Af-|- l/2)S-l-?i/ol  is  bounded  for  all  values  of  2/„;  that  is,  there  is  a  constant  W  such  that 
l/2)S-ri;/o]|<Cll^-  The  absolute  value  of  the  integral  is  less  than 

dyo 


25(A^-1- 1/2)  IT  J_^  |(iy___.)2_!'^^V+ 1/2)26’'!  ^ 2S{A^-f  1 /2)  B  (y^_y)2+(^V-M/2)2S^-r^ 

25(A'-f  i/2)iF  _r.„„-i  _ _ til _ -.-inn-  ~i 

~  ^(N+li2yS-^-x±  y  (N+ll2yS-^-x^  V(iV+ l/2)^S^-.rJ 

By  the  residue  theorem, 

r  dZo=2rriw' (z) 
J  A  ^0  * 

The  periodicity  of  w'iz)  implies  that 


As  0=  and  yg-^d,  this  quantity  approaches  zero.  It  has 
been  shown,  therefore,  that  when  <-»“  and  ■yg—’d, 

f  ^'’^dz,-^iTi(w,'  +  w,')  (B4) 

J R  ^0  ^ 

..4X..SS  iTSb'"- 

- /..liS  ■''•+S.f..  ctSS  *-J..  s  s  <'*• 


(B5) 


(B6) 


as  A'—*  “ . 

IVhen  eipiations  (B4),  (Bo),  and  (B6)  are  substituted  into 
eiiuation  (Bl).  the  e.xpression  for  the  comple.x  velocity  is 
obtained: 

(-)  =  i  A.  J^j,  ie'(2„)  cot  ^  {z„-z)  dz,  (B7) 


The  comple.x  potential  is  obtained  from  equation  (B7)  by 
integrating  with  respect  to  s  and  neglecting  the  arbitrary 
constant, 

w{z)  =  z'wj  +  —.  [  w'(z„)  log  sin  5  (z— 2o)  dzo  (B8) 

-TTi  Jfio  ^ 

where  ■w„'=^'  't-—  is  the  mean  stream  velocity. 


APPENDIX  C 


SYMBOLS 


Tlu'  priiu'ipul  symbols  us<'d  tliroii^hout  the  report  arc 
listed  here  for  eoiiveiiieiiee  of  refereiiee. 

.ii=  4^  log  siir^,  (A--j-„)-|-sinh-^,  (y-y„)J 

K  natural  logarillmi  of  <listnnee  from  singular  |)oiut  to 
center  of  eirelc  eorrespoiuling  to  cascade  airfoil 

1  total  arc  length  of  airfoil 

S  distance  between  successive  airfoils  in  cascade 
s  arc-length  parameter  corresponding  to  s 
arc-length  parameter  corresponding  to  Zo 
T’  magnitude  of  uniform  or  mean  stream  velocity  in  air¬ 
foil  or  cascade  plane  (fig.  1) 

Fj  magnitude  of  uniform  stream  velocity  in  circle  ))lane 
Vi  x-component  of  uniform  or  mean  stream  velocity  T' 
Ft,  r  resultant  local  mean  stream  x-component  of  vidocity  F 
Fj  )/-component  of  uniform  or  mean  stream  velocity  F 
Vc  local  v('locit3'  on  circle  corresponding  to  isolated 
airfoil 

Vc.c  local  velocity  on  circle  corresponding  to  airfoil  in 
cascade 

«T  velocity  induced  by  vortices  in  region /2>0 
w  complex  potential  function, 

wj  complex  velocity  of  mean  stream  for  airfoil  in  cascade 

W+w,')  =  Vi-iV^ 

u\/  complex  velocitv*  of  uniform-stream  for  isolated  airfoil, 
f  y 

X  real  part  of  s 

Xr.i/c  coordinates  of  point  about  which  airfoil  is  rotated 
(centroid  of  vortex  distribution  for  cascade  airfoils) 
y  imaginary  part  of  z 

2  coordinate  of  point  where  stream  function  is  com- 

putccl,  x-j-iy 

2„  coordinate  of  point  where  vortex  is  located,  Xo+iy„ 

a  angle  of  inclination  of  iiniform  stream  velocity  to 

x-axis 

/3  angle  through  which  airfoil  is  rotated 

r  circulation  about  airfoil 

y(z„)  vortex  strength  per  unit  arc  length  at  2„ 

6  central  angle  of  circle 

fl.v  angle  of  stagnation  point  on  circle  corresponding  to 
loading  edge  of  airfoil 

8t  ■  angle  of  stagnation  point  on  circle  corresponding  to 
trailing  edge  of  airfoil 


X  luigle  of  inclination  of  mean  How  to  normal  to  ca 
axis  (fig.  1) 

(fl  velocitv  potential  on  airfoil.  7f[ie(2')l 

ifi.uipii,  values  of  tp  at  points  ^1,  B,  C,  D,  respectivelv, 
<Pc,Vn  curve  of  intersects /2(.s,  x„)=0  (See  iig.  ■ 
<Pc  velocitv  potential  on  circle  corresponding  to  is( 
airfoil 

Pc,c  velocity  potential  on  circle  corresponding  to  air; 
cascade 

<Pmin  velocity  iiotential  at  leading  edge  of  airfoil 
i/'  stream  function,  I[w{z)\ 

\p,„  stream  function  of  mean  stream  of  cascade  flow 
\pu  stream  (unction  of  uniforrn  stream  flowing  . 
isolated  airfoil 

mean  value  of  stream  function  over  airfoil 
Ai/-  variation  of  stream  function,  if- 
Subscripts  1  and  2  when  appended  to  w' ,  F,  and  T'Tiiu 
inflow  and  discharge  values,  respeetively. 
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TABLE  1.  C’OORDIXATKS  OE  Ofx  -.r„.  // 
(a)  Valiu's  {)f  (!f  —  ij„)/S 


(1 


0. 0*1 


-(1.40  i  o.n2:>7 

-.:{S  1  .01102 

i  .o;{;u 

;  .o:c:> 

-,;52  i  .012*1 


0.  10 


I 

I 


0.  15 


I 

0. ;«)  1 


0.  10  (1,  4‘)  0.  •)(l 


-0.  ;m  i  0  (MSI 


-.•28 

1  . 0,545 

.  0.0229 

1 

! 

-  20 

1  .0018 

.0371 

1 

-.24 

j  . 0009 

.  0497 

1 

1 

-.22 

1  .0791 

.0021 

-0.  '20 

0.  OS'J  I 

0.  07.50 

1 

j 

-.  18 

.  1010 

.IWHT 

1  n,9-290 

1 

i 

! 

1 

10 

.  1140 

1  . 1035 

1  , 0020 

1 

I 

! 

14 

.  I2S0 

i  .  1 1'.« 

1  .O.vl 

1 

j 

-.  12 

.  1447 

t  . 1309 

1  .1107 

0.0392 

j 

1 

1 

-0.  10 

0. 1020 

0. 1500 

!  0.  1344 

0. 0S81 

i 

-.08 

.  1S24 

.  iras 

1  ,  l.'.8K 

1241 

0.0454 

! 

1 

- .  00 

.  2011 

.  1993 

1  ,  1.844 

.  1572 

.  103 

-.04 

'.>■>77 

.  2-230 

:  .2113 

.  18% 

.  1558 

0.  1014 

{ 

1 

-.02 

.  2532 

.  2IU8 

i  . 2390 

.  2222 

.  IIMHI 

.  1608 

0. 1090 

1 

0 

0.  280.5 

0. 2778 

0-  2094 

0,  2.«.1 

0.  23.54 

0. 2(HH) 

0. 1777 

0.  14(H) 

1).  0909 

(1. 1)190 

1  0.  ()(KH) 

.02 

.  3097 

.  3074 

.  3000 

-  2.892 

.  27:i7 

.  2542 

.  2318 

,aisi 

.  1858 

.  1092 

1  . 1029 

.04 

.  340.5 

.  3380 

.  3331 

-  3239 

.3117 

.  *2908 

.  2.^04 

.  2038 

,2191 

,  2389 

1  . 23.52 

.00 

.3728 

.3713 

.  3r4>K 

.  3.595 

.3498 

.  3:484 

.  :4200 

.3139 

.  3031) 

.  2!H’)5 

.  2940 

.08 

.  4U(M 

.  40.52 

.4010 

.  39.58 

.3881 

.3793 

.  :4698 

,  31)08 

.  3.533 

.  :i482 

.:}164 

0. 10 

0. 4412 

0. 4402 

0. 4373 

0.  4327 

0. 4207 

O.fllOS 

0.4120 

0,  4058 

0,4001 

0, 3904 

■0. 3951 

.  12 

.  4709 

.  4701 

.  it:i9 

-  4702 

.  40.55 

..4001 

.  4,545 

,  149;4 

.4451 

.  14*23 

.441:4  i 

.  14 

..5135 

.  .5r29 

.5111 

.  5082 

.  5(M5 

.■5003 

.  1900 

.  49*20 

,  4888 

.  4807 

.  48.59 

.  10 

.  ,5.508 

.  .5503 

.5488 

.  5460 

.  5437 

.'.5404 

.  5371 

.  5:440 

,,5:410 

.  5299 

.  .5294  1 

•  .  18 

.  5880 

.  .5882 

-  .5.871 

..5853 

.5830 

.  .5805 

.  5779 

.  57jJ5 

.  5730 

.  5724 

,  57*20 

0,  20 

0, 0209 

0,  0200 

0.  0257 

0. 6243 

0.6*225 

0.6205 

0.  0185 

0, 0107 

0.  0152 

0. 0143 

0,  0140 

.2'2 

,  005.5 

.  00.53 

.  CiiMO 

.  66:45 

.  6621 

.  (HMO 

.  0.590 

.  0576 

.  0505 

.  0557 

.  (’>555 

.24 

.  7044 

.7012 

.  7037 

.7029 

.7018 

.7006 

.  0!W4 

.  t'.983 

.  0974 

.  0(H)8 

.  0900 

.  'iO 

.  7430 

.7434 

.7430 

.7424 

.7415 

.  74(H) 

.  7:490 

.  :.'88 

.7:481 

“37" 

.  7375 

.28 

.7829 

.7828 

.  7825 

.7820 

.7813 

.7800 

.7798 

.7792 

.7787 

.7783 

.  7782 

0.  30 

0.  8224 

0.  8223 

0. 8221 

0.  8217 

0. 8211 

0. 8206 

0,  8*200 

0.8195 

0,8191 

0.8188 

0.8187 

,32 

.  Sti'iO 

.  8019 

.8017 

.8614 

.  8610 

.8606 

.8001 

.  8.597 

.  8,594 

.  8592 

.  8592 

.34 

.9017 

.  9010 

.9015 

.9012 

.9009 

.9000 

,9002 

.  8999 

.  8997 

.8996 

.8995  1 

.3(1 

.9415 

.9414 

.9413 

.9411 

.9409 

.  9400 

.9403 

.9401 

.  9:499 

.  9398 

. 9397  j 

.38 

.9813 

.  9812 

.9811 

.9810 

.9808 

.0800 

.  9804 

.  9802 

,  9800 

.9800 

.9799  1 

0.40 

1.0211 

1.0211 

1.0210 

1.0200 

1.0207 

1.0206 

1.0204 

1 , 0293 

1.0*201 

1.0201 

1.0201 

.42 

1.0010 

1.0009 

l.(Hi08 

l.(Mi07 

1.4M400 

\.m)A 

l.or>03 

1.000*2 

1 . 0i)02 

I . 0002  1 

.44 

1.  l(K)U 

1.  10(19 

1.  KNIH 

1.  KXW 

1.1007 

i.lO(M4 

1.  1005 

.  1.  1(X)4 

1,  HH):4 

1.  1(X)3 

1.  1003  1 

40 

1.  1408 

1.  14(18 

1.  1408 

1.  1407 

1.1407 

I.HIMI 

1.  1405 

1.  1404 

1, 1404 

1.  140:4 

1. 1403 

.48 

1.  1808 

1,  1807 

1.1507 

1.  1807 

1.  1806 

1.1800 

1.  1805 

1.  1805 

1.  1804 

1,  18(H 

1.  1804 

0.50 

1.2207 

1,2207 

1.2207 

1.  22(Hi 

1. 2*206  1 

1.  2*200 

1.  2*205 

1,2*205 

1,2204 

1  2'204 

1.2*204 

52 

1.2C>07 

1.2007 

1. 2007 

1.  iK'itHi 

1.2606  1 

.  1.26(H4 

1.  2r)()5 

1,2005 

l.*20()5 

i.  2005 

1,2005 

..'■>4 

1 . 3007 

1.3000 

1.3000 

1.3006 

I.30U6  1 

1.  :4(NH'> 

l.;4005 

1,  :41H)5 

1.3(K)5 

l.;4(X)5 

1,:4005 

.  50 

1.3400  1 

1.3400 

1.3406 

1.34% 

l.:4406  1 

1.34% 

l.:4405 

1.;J405 

1.3405 

l.;4405 

l.:44()5 

..58 

1.3S0G  1 

1.3800 

l-ltSOC  ! 

i.:4sor> 

1.3806 

i.:4800  1 

1.  :4S05 

l.:4,805 

1 , 3805 

l.;5Sl)5 

1.3805 

1).  00  1 

1.4200  i 

1.4200 

1.4206  1 

1.42(H)  i 

1.4*206  . 

1.4200  1 

1.4*205  1 

1.  4205 

1,4205 

1 ,  4*205 

1.4205 

0-2  ! 

1.4t'i00  ! 

l.KWO 

I  1006 

1 .  40(H>  1 

l.KHH)  i 

1. 4(>00  1 

l.4%5  ! 

1.  4(H)5 

1 ,  0H)5 

1 ,  4(H)5 

1,  4)'r()5 

04 

1.5(K10  , 

1.5iH)0 

i..*8xm 

1 .  .500'i  ! 

1.5(H)6  < 

1.5000  1 

1..5(Hk5 

1,  .5(M)5 

i .  ,•005 

1,  5(K)5 

1.  .5005 

00  1 

1..54(Hi 

1..5H)0 

1..54(H) 

1.5404)  1 

1.54(M)  j 

l.,54(«l  1 

1.  5400 

1 . 5405 

1..540,5 

1,  5105 

1 . 5405 

.08  [ 

1.3800  r 

1,5800  . 

1 . 5806 

1.5806  ! 

1.58(Hi  j 

1.5800  i 

1.5800  ' 

1.OS05 

1.  5805 

1.58U5 

1.  .5805 

0.  70 

1.0'200 

1.0206  , 

1.0206 

I.020C  1 

1 

1.0206  j 

1.6200  1 

1 

1.6*200 

1.0205 

1,  0205 

1.0205 

1.0205 

TABLE  I.  COORDIXATES  OF /sU-x,,  //-y,)— Concluded 

(1))  Values  of  (x— xJ/S 


!\  v-p.  1  '  j 

■  \  '“s’  !  I  j 

1  \ - >  0  0.025  U.0.M)  0.075 

i  "!  \  1  ! 

0. 100 

0.  125 

0.  150 

'  '  i 

0.175  0.  2(K)  0.*225  j  0.250 
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i 
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.0782 
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0.1137  0.0685 
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.2317  .'MiH 
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APPENDIX  C 


SYMBOLS 


Tlio  priiu'ipal  symbols  usml  tbroiiirhout  tlio  l■(■po^•t  are 
listed  here  for  eouvenieiiee  of  reference. 

ji=  (-f  -  J-J J 

log  siir  -(-sinlr  (i/— i/„)J 

K  natural  logarilluu  of  distance  from  singular  point  to 
center  of  circle  corres()onding  to  cascade  airfoil 
I  total  arc  Icirgth  of  airfoil 
(S'  distance  between  successive  airfoils  in  cascade 
s  arc-length  parameter  corresponding  to  r 

Sa  arc-length  parameter  corresponding  to  z„ 

F  magnitude  of  uniform  or  mean  stream  velocity  in  air¬ 

foil  or  cascade  plane  (fig.  1) 

Fj  magnitude  of  uniform  stream  velocity  in  circle  plane 
Fi  j:-compoiient  of  uniform  or  moan  stream  velocity  F 
Fi,  r  resultant  local  mean  stream  a:-component  of  vedoeity  T' 
F„  y-component  of  uniform  or  mean  stream  velocity  F 
Vc  local  velocity  on  circle  corresponding  to  isolated 
airfoil 

t’e,  c  local  velocit}^  on  circle  corresponding  to  airfoil  in 
cascade 

w,  velocity  induced  by  vortices  in  region 

w  comple.v  potential  function,  (p+itp 

w„'  comple.x  velocity  of  mean  stream  for  airfoil  in  cascade 

(w/+te2')  =  Fx— iFjJ 

comple.x  velocity  of  unifornrstream  for  isolated  airfoil, 

T  X  1 1 1/ 

X  real  part  of  z 

coordinates  of  point  about  which  airfoil  is  rotated 
(centroid  of  vortex  distribution  for  cascade  airfoils) 
y  imaginary  part  of  s 

z  coordinate  of  point  where  stream  function  is  com¬ 
puted,  x-\-iy 

z„  coordinate  of  point  where  vortex  is  located,  Xo+bjo 
a  angle  of  inclination  of  uniform  stream  velocity  to 
£-axis 

/3  angle  through  which  airfoil  is  rotateil 

r  circulation  about  airfoil 

7(2„)  vortex  strength  per  unit  arc  length  at 
6  central  angle  of  circle 

8,v  angle  of  stagnation  point  on  circle  corresponding  to 
leading  edge  of  airfoil 

6t  ■  angle  of  stagnation  point  on  circle  corresponding  to 
trailing  edge  of  airfoil 


X  angle  of  inclination  of  mean  How  to  normal  to  ca 
axis  (fig.  1) 

<p  velocity  potential  on  airfoil.  /f[ic(r')] 

V.i.Vs,  values  of  (p  at  points  ^1,  B,  C,  I),  respectively, 
ipc,<PD  curve  of  <p(s„)  intersects /^(s,  .«„)=0  (See  fig. 
ipe  velocity  potential  on  circle  corresponding  to  is( 
airfoil 

ipe.c  velocity  potential  on  circle  corresponding  to  air; 
cascade 

fmin  velocity  potential  at  leading  eilge  of  airfoil 
\p  stream  function,  I[w{z)] 

tpm  stream  function  of  mean  stream  of  cascade  flow 
stream  function  of  uniform  stream  flowing  ; 
isolated  airfoil 

mean  value  of  stream  function  over  airfoil 
variatfon  of  stream  function,  ^ 

Subscripts  1  and  2  when  appended  to  w' ,  F,  and  F,  im 
inflow  and  discharge  values,  respeetively. 
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ISOLATED  A.\D  CASCADE  AIHFOILS  WITH  PUESCKIBED  VELOCITY  DISTRIBUTIOX 
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TAIJIJ':  ].  t'OOlMHXATES  OF /,(x  -.r„. 
(;li  Values  of  (i/  — ;/„)/S 
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i  ; 

1  i 
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0.  1020 

0.  1.500 

1  0.  1344 

0.  OSSl 

i 

j 

1 

-.OS 

.  1.824 

,  170.8 

.  1.5SH 

.  1211 

0. 04.54 

1 

i 

1 

1 

-.oil 

.20)1 

.  1993 

.  1.844 

.  1.572 

.1103 

! 

-.04 

'.'•.>77 

.  2230 

.2113 

.  1S90 

.  1.558 

0. 1014 

! 

1 

-.02 

.  2532 

.2198 

.  2390 

.  2222 

.  1900 

.  IIKJH 

0. 1090 

( 

0 

0. 2805 

0,  27Ta 

0,  2094 

0.  2.553 

0.  23.54 

0.  2090 

0.  1777 

0,  141M) 

0. 0909 

0.0490 

1  0, 0000 

.02 

.  3007 

.  3074 

.  3000 

.  2.892 

.  2737 

.  2.542 

.2318 

.2081 

.  18.58 

.  1092 

1  . 1029 

.04 

.  3405 

.  3380 

.  3331 

.  3239 

.3117 

.  2908 

.  2.^04 

.  2038 

,  2-lUl 

,  2389 

i  . 2352 

.0(1 

.3728 

.3713 

.  30)W 

.  3.595 

.3498 

.  3384 

.  3200 

.3139 

,  3030 

,  2905 

.  2940 

.08 

.  4004 

.  40.52 

.  4010 

.  39.58 

.3881 

.3793 

.  3098 

,  3008 

.  3.533 

.  3482 

.3104 

0. 10 

0.  4412 

0.4402 

0. 4373 

0.  4327 

0. 4207 

0.)ll9« 

0.4120 

0. 4058 

0,4001 

1),  39tM 

0.3951 

.  12 

.  1701) 

.  4701 

.  I7:iu 

.  4702 

.  4055 

..4001 

.  4.545 

,  4403 

.4451 

,  1423 

.4413 

.  14 

..5135 

.  .5129 

.5111 

.  .5082 

.  ,5(M5 

■5003 

.  49»K) 

.  4920 

.4888 

.  4807 

.  48.59 

.  10 

.  .5.*0.S 

..5503 

.  .5488 

.  .5400 

.5437 

.5104 

.  5.371 

,  ,5340 

.  .5310 

.  ,5299 

.  5294  i 

•  .  18 

.  .5880 

.  .5.882 

.  .5871 

.  58.53 

.5830 

.  5805 

.  5779 

.  5730 

,  5724 

.  5720 

0.20 

0,  0209 

0.  0200 

0.  0257 

0.  0243 

0.0225 

0-  0205 

0.  0185 

0. 0107 

0,0152 

0.  0143 

0. 0140 

.22 

.  0055 

.  00.53 

.  0040 

.  0(k15 

.  0021 

.  (KiOO 

.  (1590 

,  0.570 

.  0.505 

.  0,557 

.  0555 

.24 

.  7044 

.  7012 

.  7037 

.7029 

.7018 

.  7(M)0 

.0994 

.  0983 

.  0974 

.  01H)8 

.  0900 

.20 

.  7430 

,7434 

.7430 

.7424 

.7415 

.  7400 

.  7390 

.  7.588 

.7381 

.7377 

.  7375 

■=s._ 

.7829 

,  7828 

.  7825 

.7820 

.7813 

.7800 

.7798 

.7792 

.7787 

,7783 

. 7782  1 

0.  30 

0.8224 

0.  8223 

D.  8221 

0.  8217 

0. 8211 

0.8200 

0-  8200 

0.8195 

0,8191 

0.8188 

0.8187 

.32 

.  8020 

.  8019 

.  S0l“ 

.8014 

.8010 

.8000 

..H(K)1 

.  s.5y7 

.  8594 

.8592 

.8592 

.34 

,9017 

.9010 

.9015 

.9012 

.9009 

.9000 

.  1HK)2 

.  8999 

,8997 

.8995 

.8995  , 

.30 

.9415 

.9414 

.9413 

.9411 

.9409 

.  ‘MOO 

.9403 

.  9401 

.  9399 

.  9398 

. 9397  j 

.38 

.9813 

.9812 

.98U 

.9810 

.9808 

.9800 

.9804 

.  9802 

.  ilWIO 

.9800 

.9799  ' 

0.  40 

1.0211 

1.02U 

1.0210 

1.0209 

1.0207 

1.0200 

1.0204 

1 . 0203 

1,0201 

1.0201 
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.42 
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1.  1009 

1.  1009 
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1.2005 
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1.3005 
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i.awii  1 
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1.3S06  1 

1.3800 
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l..541Hi  1 
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1.  5805 
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0.  70 

1.0200  1 

1.020G 

1.0200 
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1 
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1.0200  1 
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TAI3LK  I.  (.'OOUniXvVTKS  OK  /2(x— u^Vo) — ConcliKlod 
(b)  Values  of  {x-~Xo)!S 
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10 
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.  I33(i 
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Y 


X 


Positive  directions  of  axes  and  angles  (forces  and  moments)  are  shown  by  arrows 


Axis 

Force 
(parallel 
to  axis) 
symbol 

Moment  about  axis 

Angle 

Velocities 

Designation 

H 

Designation 

S3rm- 

bol 

Positive 

direction 

Designa¬ 

tion 

Sym¬ 

bol 

Linear 
(compo¬ 
nent  along 
axis) 

Angular 

T.nngitiiHinal 

X 

X 

Rolling . . 

L 

Y - >Z 

Roll  .  .. 

u 

P 

Lateral _  _ 

Y 

Y 

Pitching. _ 

M 

Z - tX 

Pitch 

e 

V 

q 

Normal . 

z 

Z 

YawingT _ 

N 

X - *Y 

Yaw _ 

w 

r. 

Absolute  coefficients  of  moment 

c-Jl-  O-^  c--^ 

(rolling)  (pitching)  (yawing) 


Angle  of  set  of  control  surface  (relative  to  neutral 
position),  5.  (Indicate  surface  by  proper  subscript.) 


4.  PROPELLER  SYMBOLS 


D  Diameter 

p  Geometric  pitch 

pjD  Pitch  ratio 

V'  Inflow  velocity 

V,  Slipstream  velocity 

.  T 

T  Thrust,  absolute  coefficient  (7t= — 57=3 

prrD* 

Q  Torque,  absolute  coefficient  C7g=— ^ 


P 

P 

Power,  absolute  coefficient  (7p= — 5™ 

c. 

5 

Speed-power  coefficient 

V 

Efficiency 

n 

Revolutions  per  second,  rps 

\2vrn/ 

Effective  helix  angle = tan”' 

6.  NUMERICAL  RELATIONS 


1  hp=76.04  kg-m/s=550  ft-lb/sec 
1  metric  hor8epower=0.9863  hp 
1  mph=0.4470  mps 
I  mpB= 2.2369  mph 


1  lb=0.4536  kg 
1  kg=2.2046  lb 
1  mi=  1,609.35  m=5,280  ft 
1  m=3.2808  ft 


